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Abstract 

Background: Salmonids are popular sport fishes, and as such have been subjected to widespread stocking 
throughout western North America. Historically, stocking was done with little regard for genetic variation among 
populations and has resulted in genetic mixing among species and subspecies in many areas, thus putting the 
genetic integrity of native salmonid populations at risk and creating a need to assess the genetic constitution of 
native salmonid populations. Cutthroat trout is a salmonid species with pronounced geographic structure (there are 
10 extant subspecies) and a recent history of hybridization with introduced rainbow trout in many populations. 
Genetic admixture has also occurred among cutthroat trout subspecies in areas where introductions have brought 
two or more subspecies into contact. Consequently, management agencies have increased their efforts to evaluate 
the genetic composition of cutthroat trout populations to identify populations that remain uncompromised and 
manage them accordingly, but additional genetic markers are needed to do so effectively. Here we used genome 
reduction, MID-barcoding, and 454-pyrosequencing to discover single nucleotide polymorphisms that differentiate 
cutthroat trout subspecies and can be used as a rapid, cost-effective method to characterize the genetic 
composition of cutthroat trout populations. 

Results: Thirty cutthroat and six rainbow trout individuals were subjected to genome reduction and 
next-generation sequencing. A total of 1,499,670 reads averaging 379 base pairs in length were generated by 454- 
pyrosequencing, resulting in 569,060,077 total base pairs sequenced. A total of 43,558 putative SNPs were identified, 
and of those, 125 SNP primers were developed that successfully amplified 96 cutthroat trout and rainbow trout 
individuals. These SNP loci were able to differentiate most cutthroat trout subspecies using distance methods and 
Structure analyses. 

Conclusions: Genomic and bioinformatic protocols were successfully implemented to identify 125 nuclear SNPs 
that are capable of differentiating most subspecies of cutthroat trout from one another. The ability to use this suite 
of SNPs to identify individuals of unknown genetic background to subspecies can be a valuable tool for 
management agencies in their efforts to evaluate the genetic structure of cutthroat trout populations prior to 
constructing and implementing conservation plans. 

Keywords: Conservation genetics. Genetic admixture. Hybridization, KASPar, Oncorhynchus clarkii, Oncorhynchus 
mykiss, Population genomics, SNP, Rainbow trout 
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Background 

Single nucleotide polymorphisms (SNPs) are powerful 
genetic markers that are increasingly being used in phylo- 
genetic and population genetic studies [1,2]. Advances in 
high-throughput sequencing technologies have made SNP 
identification faster and cheaper than traditional methods 
that utilize Sanger sequencing (see Metzker [3] for a review 
of next-generation techniques). The relative ease by which 
SNPs can now be identified makes SNP discovery more at- 
tainable [2,4]. Indeed, many researchers have recently used 
next generation sequencing to detect SNPs in a variety of 
organisms [5-9]. SNP discovery in salmonid fishes has gar- 
nered much attention in recent years [10-19], but there is a 
growing need for additional SNP discovery in certain 
groups of salmonids. 

Cutthroat trout, Oncorhynchus clarkii, a native western 
North American salmonid, has ten extant subspecies (and 
two extinct subspecies). The species appears to be mono- 
phyletic [20], and diversification among subspecies is postu- 
lated to have begun approximately two million years ago 
[20,21]. The extant subspecies are: Bonneville cutthroat 
trout (O. c. Utah) in the Bonneville Basin; Coastal cutthroat 
trout (O. c. clarki) in coastal drainages from Alaska to 
northern California; Colorado River cutthroat trout (O. c. 
pleuriticus) in the upper Colorado River basin; Greenback 
cutthroat trout (O. c. stomias) in the Arkansas and the 
South Platte river basins in eastern Colorado (federally 
listed as threatened); Lahontan cutthroat trout (O. c. 
henshawi) in the western Lahontan Basin of Nevada and 
several closed basins in Oregon (federally listed as threa- 
tened); Humboldt cutthroat trout (O. c. humholdtemis) in 
the Humboldt River in the eastern portion of the Lahontan 
Basin recently designated as a separate sub-species from 
the Lahontan cutthroat trout [22]; Paiute cutthroat trout 
(O. c. seleniris) in Silver Creek on the eastern slope of the 
Sierra Nevada (federally listed as threatened); Rio Grande 
cutthroat trout (O. c. virginalis) in tributaries to the Rio 
Grande in southern Colorado and New Mexico; Westslope 
cutthroat trout (O. c. lewisi) in drainages of the Rocky 
Mountains in Alberta, British Columbia, northern Idaho, 
and Montana, with disjunct populations in Oregon and 
Washington; Yellowstone cutthroat trout (O. c. bouvieri) in 
the Yellowstone River, Yellowstone Lake, and the upper 
Snal<e River drainages of Wyoming, Idaho, and Montana. 

The Snake River fine spotted cutthroat trout (in the 
upper Snake River) was designated as a separate subspecies, 
O. c. behnkei [23]. The only documented differences be- 
tween it and sympatric Yellowstone cutthroat trout are the 
spotting pattern and behavior, whereas other morphological 
and meristic characters are the same [23,24]. Genetic ana- 
lyses have not revealed differences between Snake River fine 
spotted and Yellowstone cutthroat trout [21,25,26]. Unfor- 
tunately, Montgomery [23] did not designate a type speci- 
men when he named O. c. behnkei, thus technically ren- 



dering the subspecies designation invalid. Therefore, we 
treat the Yellowstone cutthroat trout and the Snake River 
fine spotted cutthroat trout as O. c. bouvieri. Similarly, 
while Paiute cutthroat trout have fewer spots than 
Lahontan cutthroat trout (many have no spots), they do 
not differ firom Lahontan cutthroat trout with any other 
morphological or meristic characters that have been exam- 
ined [24], nor do the two subspecies appear to be genetic- 
ally distinct based on electrophoretic data [27]. Recently, 
Finger et al. [28] characterized six SNPs that differentiate 
rainbow trout from Lahontan and Paiute cutthroat trout, 
but the two cutthroat trout subspecies were identical at all 
six SNP loci. Additionally, Paiute cutthroat trout and 
Lahontan cutthroat trout carried identical haplotypes 
of the second subunit of the NADH dehydrogenase- 
ubiquinone oxidoreductase enzyme complex I (ND2) of 
the mitochondrial genome [21]. Hence, Paiute cutthroat 
trout and Lahontan cutthroat trout are likely to be very 
similar genetically. 

The Bonneville cutthroat trout includes a morphologic- 
ally and ecologically unique lineage of cutthroat trout in the 
Bear River drainage [24,29]. Based on allozyme [30,31] and 
mtDNA [21,26,32] data, the Bear River strain of Bonneville 
cutthroat trout (hereafter referred to as Bear River cut- 
throat trout for simplicity) are genetically more close- 
ly related to Yellowstone cutthroat trout than to the 
Bonneville cutthroat trout in the main Bonneville Basin 
[26,30-32], although an associated taxonomic revision has 
not yet been made. The sister relationship between 
Yellowstone and Bear River cutthroat trout makes biogeo- 
graphic sense because the Bear River was part of the upper 
Snake River drainage until the late Pleistocene, at which 
time the Bear River was redirected into the Bonneville 
Basin (~35 Ka) [33-35]. Contemporary gene flow between 
the Bear River and other drainages within the Bonneville 
Basin is prevented by the Great Salt Lake. While the cut- 
throat trout in the Bear River system are currently classified 
as Bonneville cutthroat trout, herein we treat them as an 
additional lineage within the species O. clarkii. 

Despite these unique lineages, cutthroat trout were 
stocked within and among major drainages with little con- 
cern for genetic variability among subspecies [24,36]. Intro- 
gression among cutthroat trout subspecies has resulted 
from these stocking practices. Additionally, rainbow trout 
{Oncorhynchus mykiss) have been stocked extensively 
throughout western North America. Rainbow trout readily 
hybridize with cutthroat trout in areas where the two spe- 
cies did not formerly co-occur, posing a serious threat to 
the genetic integrity of the native cutthroat trout popula- 
tions [25,37-39]. Rainbow trout x cutthroat trout hybrids 
can be identified reasonably well using morphological and 
meristic characters, however in populations with extensive 
introgressive hybridization this is not always the case, nor 
are hybrids between cutthroat trout subspecies easily 
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recognized. Because of widespread intra- and interspecific 
hybridization, management agencies have increased efforts 
to assess the genetic composition of native cutthroat trout 
populations. A range of genetic markers have been used to 
assess introgressive hybridization in cutthroat trout popu- 
lations, and recently SNPs useful for species identification 
and detecting introgression between rainbow trout and 
cutthroat trout have been developed [28,40,41]. Some 
SNPs differentiate certain subspecies of cutthroat trout 
from rainbow trout (e.g., rainbow trout vs. westslope cut- 
throat trout [39,42]; rainbow trout vs. westslope and 
Yellowstone cutthroat trout [43]; rainbow trout vs. west- 
slope, Yellowstone, coastal, and Lahontan cutthroat trout 
[44]; rainbow trout vs. westslope, Bonneville, Yellowstone, 
coastal and Lahontan cutthroat trout [45]). However, it 
remains unclear whether these SNPs are unique to those 
cutthroat trout subspecies, or if they are shared with other 
subspecies not included in those studies. Here we use 
next-generation sequencing technologies to identify add- 
itional nuclear SNPs that collectively characterize most of 
the cutthroat trout subspecies. Specifically, we used gen- 
ome reduction, MID-barcoding, and 454-pyrosequencing 
in an attempt to discover nuclear SNPs that differentiate 
nine lineages of cutthroat trout (i.e., Bear River cutthroat, 
Bonneville cutthroat, coastal cutthroat, Colorado River cut- 
throat, greenback cutthroat, Lahontan cutthroat, Rio 
Grande cutthroat, westslope cutthroat, and Yellowstone 
cutthroat) and rainbow trout (including Columbia redband 
trout and steelhead). These SNPs were used to develop a 
SNP assay that can be easily used to evaluate the genetic 
integrity of cutthroat trout populations across the entire 
range of the species. The SNP assays are based on 
KASPar genotyping chemistry and were detected using 
the Fluidigm dynamic array platform. 

Methods 

A brief overview of our methods is as follows: DNA was 
extracted from tissue samples, followed by genome reduc- 
tion and 454 pyro-sequencing. The sequences produced 
were assembled and scanned for SNPs bioinformatically. 
SNPs were then genotyped and SNP diversity analyses 
were performed. Individuals that yielded unexpected 
results according to a priori subspecies designations were 
re-sequenced for the ND2 mitochondrial DNA gene using 
Sanger sequencing for further evaluation. 

DNA extraction 

Fin clips or muscle tissues were obtained from thirty-six 
individuals representing nine sub-species of cutthroat trout, 
plus the Bear River cutthroat trout and rainbow trout 
(including Columbia redband trout and steelhead, which 
are unique forms of O. mykiss) from the Monte L. Bean Life 
Science Museum ichthyological collection (Table 1). The 
samples were collected by field biologists familiar with 



cutthroat trout subspecies identifications based on pheno- 
typic characters. While subspecies identifications are typic- 
ally accurate when made by experts, especially when 
accounting for geographic distribution, the presence of 
cryptic hybrids among these samples is possible. We were 
unable to obtain tissues for Paiute cutthroat trout, but given 
the lack of distinguishing genetic characteristics between 
them and Lahontan cutthroat trout, they are likely very 
similar to Lahontan cutthroat trout. Whole genomic DNA 
was extracted using Qiagen DNeasy Tissue kits following 
the manufacturers recommended protocol. All extracted 
DNA was quantified using a NanoDrop 1000 Spectropho- 
tometer (NanoDrop Technologies Inc., Montchanin, DE), 
and each sample was diluted to a concentration of 150 
ng/|aJ using nuclease free water. 

Some of the DNA samples that were extracted were not 
high enough quality to be included in the genomic reduc- 
tion and 454-pyrosequencing steps (outiined below). In 
most cases this resulted in simply using one less individual 
from a subspecies, but was more problematic for Lahontan 
and Humboldt cutthroat trout where only one individual 
of each subspecies had DNA of a high enough quality to 
include. To ensure getting enough reads, those two indivi- 
duals were pooled, resulting in a Lahontan Basin complex 
where Lahontan and Humboldt cutthroat were treated as a 
single unit. Hereafter we refer to both groups as Lahontan 
Basin cutthroat, but this is merely for convenience, and is 
not meant to be a statement regarding a taxonomic update. 

To verify that the discovered SNPs were able to differen- 
tiate subspecies of cutthroat trout from each other and 
from rainbow trout, DNA was extracted from an additional 
sixty cutthroat trout individuals representing ten trout 
lineages. These sixty samples, along with the thirty-six 
samples that were extracted for 454-pyrosequencing, were 
used to verify SNPs. Hence, each cutthroat trout lineage 
was represented by ten individuals, and rainbow trout was 
represented by six individuals for SNP verification. In an 
attempt to account for geographic variation within cut- 
throat trout subspecies, samples were included from popu- 
lations different from those included on the initial SNP 
discovery panel when possible (Table 1). Samples were 
chosen from populations that are believed to be non- 
admixed, but this was not verifiable in all cases. 

Genomic reduction and 454-pyro5equencing 

We followed the genomic reduction methodology de- 
scribed in detail by Maughan et al. [6]. In brief, genomic 
reduction was carried out using restriction enzymes £coRI 
and Bfal to double digest genomic DNA at restriction sites 
that were conserved across all sub-species, and then 
attaching £coRI and Bfal adapter sequences to the sticky 
ends using T4 DNA ligase. Approximately 90% of the 
genome was then discarded through size exclusion via 
spin chromatography and biotin-streptavidin paramagnetic 
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Table 1 Sampling Information 



Species/subspecies 


n 


BYU # 


Location 


Drainage (Basin) 


Latitude/Longitude 


Cutthroat Trout (Oncorhynchus 
darkii) 












Bonneville cutthroat trout 


10 


33807 


Chalk Creek, UT 


Weber River (Bonneville) 


40° 59' 44" N, 1 1 1° 03' 31" 
W 


0. c. Utah 




134304 


Trout Creek, UT 


Deep Creek Mts. (Bonneville) 


39° 45' 55" N, 113° 55' 45" 
W 






135225 


Birch Creek, UT 


Beaver River (Bonneville) 


38° 12' 59" N, 112° 31' 49" 

VV 






1 79753 


NGDO LreeK, u i 


Spanish Fork River 
(Bonneville) 


3l U4 N, 1 11 6o ZD 

W 






93975 


Main Creek, UT 


Provo River (Bonnevi le) 


/in" 1 o' m" M 111° in' r\(^" \M 
4U ly Uz N, 1 1 1 zU Do Vv 






135133 


Willard Creek, NV 


Spring Valley 


39° or 22" N, 114° 23' 33" W 






134577 


Lost Creek, UT 


Sevier River (Bonneville) 


38° 43' 14" N, lir42' 06" W 






1 81 735 


Deadman Creek, NV 


Smith Creek (Bonneville) 


39 1 9 40 N, 1 1 4 Oo z3 W 






181774 


Pine Creek, NV 


Spring Valley 


38° 59' 09" N, 114° 22' 34" W 






181752 


Hampton Creek, NV 


Hampton Creek (Bonneville) 


39° 14' 39" N, 114° 06' 28" W 


Bear River Bonneville cutthroat 
trout 


10 


135526 


Yellow Creek, UT 


Bear River 


40° 58' 58" N, 110° 57' 35" 
W 


0. c. Utah 




136684 


Dog Creek, WY 


Bear River 


43° 37' 52" N, 110° 48' 11" 
W 






136733 


Cottonwood Creek, WY 


Upper Salt River (Bear R.) 


42° 37' 52" N, 110° 50' 34" 

W 






1 7963 1 


bmitn s rorK, wr 


Bear River 


4z lU o4 N, 1 lU DD 2.2. 
W 






135528 


Yellow Creek, UT 


Bear River 


40° 58' 58" N, 110° 57' 35" W 






1 79602 


uaniei risn natcnery, VVY 


M /A 

N/A 








\ /yooz 


Raymond Creek, WY 


Bear River (Bonneville) 


4z lo 3Z IM, 1 1 1 UU 33 VV 






1 79702 


Woodruff Reservoir, UT 


Bear River (Bonneville) 


A'\° ~)~!' CI" M 1 1 1° in' 17" \A/ 

41 1/ 3l N, 1 1 1 1 y z/ Vv 






1 79703 


Woodruff Reservoir, UT 


Bear River (Bonneville) 


A^° ~)~!' CI" M 1 1 1° in' 17" \A/ 

41 Z/ Dz N, 1 1 1 1 y z/ vv 






ziyj 1 O 


vveDer Kiver, u I 


Bonnevi le Basin 


Ac\° /ii^' c/1" M 11 n° cn' /II" \A/ 
4U 4o d4 N, 1 lU by 41 VV 


Coastal cutthroat trout 


1 n 
U 


Q1 001 

y 1 yy i 


Slippery Lake, AK 


Kuiu Island (Pacific Ocean) 




0. c. clarki 




yz/*fz 


MDernatny v_rccK, \Jn 


Willamette River (Columbia 
R.) 


4j 1 y IM, 1 ZZ 3Z 3Z 

w 






1 30929 


Lost Creek, OR 


Coos River (Pacific Ocean) 


43° 09' 48" N, 1 23° 42' 45" W 






1 30905 


Hemlock Creek, WA 


Cowlitz River (Columbia R.) 


46° 1 5' 48" N, 1 22° 43' 45" W 






1 30906 


Hemlock Creek, WA 


Cowlitz River (Columbia R.) 


46° 1 5' 48" N, 1 22° 43' 45" W 






91953 


Hackleman Creek, OR 


Willamette River (Columbia R.) 


44° 24' 45" N, 1 22° 02' 32" W 






90938 


Margaret Lake, AK 


Revillagigedo Island (Pacific) 


55° 41' 14" N, 131° 36' 24" W 






90945 


Margaret Lake, AK 


Revillagigedo Is and (Pacific) 


55 4 1 14 N, 1 31 36 24 W 






91504 


Flynn Creek, OR 


Coastal Range (Pacific) 


44° 31' 30" N, 123° 51' 44" W 






91505 


Flynn Creek, OR 


Coastal Range (Pacific) 


44° 31' 30" N, 123° 51' 44" W 


Colorado River cutthroat trout 


10 


94736 


Little Snake River, CO 


Duchesne River (Colorado R.) 


40° 59' 15" N, 107° 58' 08" 
W 


0. c. pieuriticus 




130805 


Muddy Creek, WY 


Blacks Fork River (Colorado 
R.) 


41° 10' 20' N, 110° 39' 12" 
W 






134410 


North Elk Creek, CO 


White River (Colorado R.) 


39° 54' 54" N, 107° 40' 18" 
W 






134387 


Avintaquin Creek, UT 


Yampa River (Colorado R.) 


40° 07' 26" N, 110° 44' 25" 



W 
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Table 1 Sampling information (Continued) 







135182 


White River, UT 


Price River (Colorado R.) 


39° 55' 26" N, 


111° 03' 07" W 






134518 


Cutthroat Creek, CO 


Colorado River 


37° 06' 18" N, 


1 06° 41 ' 34" W 






134537 


Navajo River, CO 


San Juan River (Colorado R.) 


37° 11' 55" N, 


1 06° 39' 49" W 






1 34658 


Steelman Creek, CO 


Williams Fork (Colorado R.) 


39° 45' 24" N, 


105° 55' 57" W 






134688 


Cabin Creek, CO 


Ranch Creek (Colorado R.) 


39° 58' 25" N, 


1 05° 44' 30" W 






1 34748 


Fryingpan River, CO 


Roaring Fork River (Colorado R.) 


39° 22' 01" N, 


107° or 58" W 


Greenback cutthroat trout 


10 


131909 


Como Creek, CO 


Arkansas River (Missouri R.) 


38° 32' 02" N, 106° 14' 17" 
W 


0. c. stomias 




132038 


South Prong Hayden Creek, 
CO 


Arkansas River (Missouri R.) 


38° 18' 47" N, 105° 49' 25" 
W 






179073 


Lower Bear Creek, CO 


Arkansas River (Missouri R.) 


38° 48' 33" N, 104° 54' 18" 
W 






180022 


Beaver Creek, UT 


Colorado River 


38° 21' 17" N, 109° 15' 48" 
W 






1 79085 


North Taylor Creek, CO 


Arkansas River (Missouri R.) 


38° 07' 14" N, 


105° 36' 15" W 






142903 


Coon Creek, CO 


Plateau Creek (Colorado R.) 


39° 05' 49" N, 


1 08° 07' 1 6" W 






1 34692 


Cabin Creek, CO 


Ranch Creek (Colorado R.) 


39° 58' 25" N, 


1 05° 44' 30" W 






1 79022 


Severy Creek, CO 


Arkansas River (Missouri R.) 


38° 52' 33" N, 


105° 02' 53" W 






78265 


West Antelope Creek, CO 


Gunnison River (Colorado R.) 


38° 35' 12" N, 


1 06° 59' 59" W 






187526 


Carr Creek, CO 


Colorado River 


39° 35' 09" N, 


108° 32' 03" W 


Lahontan Basin cuttliroat trout 


10 


96518 


Doudy Pond, UT* 


(Pyramid Lake) (Lahontan) 


41° or 17" N, 113° 58' 12" 
W 


0. c. henshawi 




96517 


Doudy Pond, UT* 


(Pyramid Lake) (Lahontan) 


41° or 17" N, 


113° 58' 12" W 


0. c. humboldtensis 




90414 


Coyote Creek, NV 


Humboldt River (Lahontan) 


40° 58' 42" N, 116° 12' 18" 
W 






1 82066 


Baker Lake, NV 


Baker Creek (Bonneville) 


38° 57' 27" N, 


114° 18' 31" W 






95362 


Virgin Creek, NV 


Alvord Basin 


41° 37' 36" N, 


1 1 9° 09' 1 5" W 






97748 


Segunda Creek, NV 


Humboldt River (Lahontan) 


40° 31' 59" N, 


1 1 5° 29' 1 0" W 






90741 


McDermitt Creek, OR 


Ouinn River (Lahontan) 


42° 00' 20" N, 


118° 00' 55" W 






91 109 


Indian Creek, NV 


Humboldt River (Lahontan) 


40° 19' 37" N, 


1 1 6° 43' 42" W 






90414 


Coyote Creek, NV 


Humboldt River (Lahontan) 


41° 19' 42" N, 


1 1 6° 33' 02" W 






91162 


IVlahogany Creek, NV 


Summit Lake 


41° 31' 01" N, 


1 1 8° 57' 24" W 


Rio Grande cutthroat trout 


10 


57036 


Torcido Creek, CO 


Rio Grande 


37° 45' 51" N, 105° 21' 35" 
W 


0. c. virginalis 




90712 


Peralta Creek, NM 


Rio Grande 


35° 41' 33" N, 106° 28' 02" 
W 






99839 


Little Vermejo, NM 


Canadian River (Mississippi 
R.) 


36° 59' 18" N, 105° 08' 04" 
W 






99858 


Middle Fork Carnero Creek, 
CO 


Rio Grande 


37° 58' 55" N, 106° 25' 32" 
W 






90714 


Peralta Creek, NM 


Rio Grande 


35° 41' 33" N, 


1 06° 28' 02" W 






99862 


IVliddle Fork Carnero Creek, CO 


Rio Grande 


37° 58' 55" N, 


106° 25' 32" W 






99844 


Little Vermejo, NM 


Canadian River (Mississippi R.) 


36° 59' 18" N, 


1 05° 08' 04" W 






57045 


Ventero Creek, NM 


Rio Grande 


36° 59' 49" N, 


105° 27' 58" W 






99841 


Little Vermejo, NM 


Canadian River (Mississippi R.) 


36° 59' 18" N, 


1 05° 08' 04" W 






90732 


Peralta Creek, NM 


Rio Grande 


35° 41' 33" N, 


1 06° 28' 02" W 


Westslope cutthroat trout 


10 


239020 


Black Bear, MT 


South Fork Flathead River 


47° 44' 41" N, 113° 23' 00" 



w 
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Table 1 Sampling information (Continued) 



0. c. Icwisi 


239024 


White Gulch, MT 


Missouri River 


46° 37' 10" N, 111° 28' 51" 
W 




187969 


Bull River, BC 


Kootenay River (Columbia R.) 


49° 33' 58" N, 115° 19' 26" 
W 




181526 


John Day River, OR 


Columbia River 


44° 16' 57" N, 118° 32' 43" 
W 




239021 


Blacl< Bear, MT 


South Fork Flathead River 


47° 44' 41 " N, 1 1 3° 23' 00" W 




239022 


IVlcGuire, IV1T 


Kootenai River 


48° 41' 20" N, 115° 16' 22" W 




239023 


IVlcGuire, IVIT 


Kootenai River 


48° 41' 20" N, 115° 16' 22" W 




239025 


WInite Gulcli, MT 


Missouri River 


46° 37' 10" N, 111° 28' 51" W 




244201 


Crazy Creek, BC 


South Thompson River 


51° 01' 43" N, 124° 33' 31" W 




244202 


Bow River, BC 


South Saskatchewan River 


51° 36' 37" N, 116° 07' 56" W 


YGllowstoRG cutthrost trout 


10 133863 


East Tsnslssp CrGGk, WY 


Bighorn River (Missouri R.) 


44° 12' 11" N, 107° 09' 53" 
W 


0. c. bouvieri 


136780 


Glade Creek, WY 


Upper Snake River 
(Columbia) 


44° 05' 40" N, 110° 43' 59" 
W 




136834 


Willow Creek, WY 


Lower Salt River 


42° 50' 45" N, 110° 52' 54" 
W 




1 36668 


Nortli Fork Horse Creel<, WY 


Upper Snake River (Columbia) 


43° 22' 51" N, 110° 37' 43" W 




141843 


Forest Creel<, WY 


Upper Snake River (Columbia) 


44° 10' 14" N, 110° 35' 49" W 




141946 


Sickle Creek, WY 


Upper Snake River (Columbia) 


44° 1 1 ' 28" N, 1 1 0° 23' 05" W 




142366 


Elk Antler Creek, WY 


Ye owstone River (Missouri R.) 


44° 36' 01 " N, 1 1 0° 28' 39" W 




90212 


Hatchery Creek, WY 


Yellowstone River (Missouri R.) 


44° 32' 57" N, 110° 24' 03" W 




179112 


Mill Creek, WY 


Yellowstone River (Missouri R.) 


44° 26' 59" N, 1 07° 27' 09" W 




179114 


Mill Creek, WY 


Yellowstone River (Missouri R.) 


44° 26' 59" N, 1 07° 27' 09" W 


Rainbow Trout {Oncorhynchus 
uiykiss) 










Rainbow Trout 


1 187968 


Rio Fuerte, Mexico 


Gulf of California 


Unknown 


Columbia Redband Trout 


3 142374 


McCoy Creek, OR 


Harney-Malheur Basin 


42° 58' 33" N, 118° 42' 32" 
W 




98251 


West Fork Bruneau River, NV 


Snake River (Columbia R.) 


41° 56' 02" N, 11 5° 40' 29" 
W 




94217 


West Fork Jarbridge River, NV 


Snake River (Columbia R.) 


41° 48' 57" N, 115° 24' 38" 
W 


Steelhead 


2 93748 


Kootenay River, BC 


Columbia River 


49° 19' N, 117° 39' W 




97060 


Midway Hatchery, LIT 


N/A 





Cutthroat trout individuals used in this study are listed along with sampling locality information for each. Individuals that were included on the 454 plate are 
highlighted with bold font. All individuals were used to validate/verify SNPs using Fluidigm nanofluidic chip arrays. 

*Transplanted population - source population was Morrison Creek, NV, which is believed to be a Pyramid Lake strain of Lahontan cutthroat. 



bead separation. Prior to 454-pyrosequencing, unique 
MID-barcodes were added to the remaining restriction 
fragments for each individual sample via polymerase chain 
reaction (PGR). DNA concentrations of the final PGR pro- 
ducts were quantified using picogreen fluorescent dye, and 
then aU PGR products were pooled in equimolar amounts 
to obtain a single sample totaling 5 ^g of DNA. Gel elec- 
trophoresis was performed on the pooled sample in a sin- 
gle lane of a 1.5% agarose gel, and DNA fragments ranging 
from -450-600 base pairs (bp) in size were removed and 



then used for 454-pyrosequencing. Sequencing was per- 
formed with Titanium reagents without DNA fragmenta- 
tion on a 454 Life Sciences Genome Sequencer FLX 
located in the DNA Sequencing Genter at Brigham Young 
University. 

Assembly and SNP discovery 

FoUowing 454-pyrosequencing, GLGBio Workbench 
bioinformatic software v.3.6.1 (Katrinebjerg, Aarhus N, 
Denmark) was used to separate the DNA sequences 
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according to the MID-barcodes that were attached to 
DNA fragments from each of the individual samples. All 
reads from all subspecies were pooled and a de novo as- 
sembly was created using Newbler v.2.6 (454 Life 
Sciences 2006-2011), after which all reads from all sub- 
species were mapped onto this assembly using the refer- 
ence mapping function in CLCBio Workbench. 

Putative autapomorphic SNPs were identified by com- 
paring sequences derived from a single cutthroat trout sub- 
species to the sequences of the other subspecies combined. 
Comparisons were completed using SNP Finder Plus, a 
custom perl script described by Maughan et al. [6]. A simi- 
lar comparison between rainbow trout and cutthroat trout 
(all subspecies sequences pooled) was used to identify spe- 
cies specific SNPs. A subset of SNPs was selected for valid- 
ation from the pool of all putative SNPs using the 
following criteria: 1) Polymorphisms were unique to a sin- 
gle subspecies of cutthroat or rainbow trout so that only 
potentially diagnostic SNPs were selected. 2) All reads at 
the SNP locus had a minimum read coverage depth > 8 to 
exclude putative SNPs that were based on too few reads. 3) 
A minimum of 50 bp existed on either side of the putative 
SNP for primer binding sites, with no indels or ambiguities 
within 20 bp of the SNP. 4) The minor allele of the SNP 
had to consist of a minimum of 3 reads comprising at least 
4% of the total alleles observed at that position or it was 
not considered as a SNP. This was an arbitrary cutoff 
designed as an extra filter to ensure that the minor allele is 
not erroneous (particularly for contigs with high coverage 
depth). 5) The SNP had at least 95% identity within the 
subspecies so that only alleles that were frxed or nearly 
fixed for a given subspecies were considered, and alleles 
that were unique to a small proportion of individuals within 
a subspecies were not. This criterion also guards against 
calling SNPs in misaligned sites. 6) The SNP did not appear 
in known repeat regions or reside within the mtDNA gen- 
ome, which was determined using RepeatMasker v.3.3.0 
[46] against the Danio reference database and by perform- 
ing BLAST searches on all contigs with a rainbow trout 
mtDNA reference genome [GenBank: NC_001717]. We ac- 
knowledge that mtDNA SNPs are useful in many cases, but 
the purpose of this study was to discover unlinked nuclear 
markers for additional statistical power in future population 
genetic studies. 7) SNPs were limited to one per contig in 
an attempt to minimize the number of linked alleles 
included. Limiting the SNPs to one per contig also served 
as a screen to eliminate paralogous sequence variants 
(PSVs). A position weighted window filter was also used to 
determine if SNPs were in a poor alignment region. Every 
polymorphism that occurred within the window was 
assigned a score based on its relative position to the SNP, 
the total score was added up and if a threshold score >6 
was reached the SNP was considered erroneous and 
skipped. This window filter was designed to avoid calling 



SNPs in misaligned regions, as well as to avoid areas that 
aligned paralogous loci, thus serving as an additional screen 
for PSVs. Primers for 288 SNP loci that met the above cri- 
teria (enough to fill three Fluidigm chips) were designed 
using the default settings in the PrimerPicker software pro- 
gram [47]. Primer sequences are listed in Additional file 1. 

SNP genotyping 

Prior to genotyping the SNPs, specific target amplification 
(STA) was used to pre-amplify each SNP locus. STA pri- 
mers are non-aUele specific and do not carry a polymorphic 
base on their 3' end (Additional file 1). STA reactions con- 
sisted of 2.5 |il of Qiagen 2X Multiplex PGR Master Mix, 
0.5 [il of the lOX STA primer assays (which consisted of 
192 [il STA primer, 192 \A constant primer [2 [il each 
primer per reaction x 96 reactions], and 16 \A TE buffer 
[10 mM Tris, 1 mM EDTA, pH 7.5]), and 0.75 ^1 DNase 
free water per reaction. A total of 3.75 \A of the STA pre- 
mix was combined with 1.25 \sl of genomic DNA and amp- 
lified with the following thermal profile: 95°G for 15 
minutes followed by 14 cycles of 95°G for 15 seconds and 
60°G for 4 minutes. Following PGR, we diluted the STA 
products (1:100) with nuclease free water. 

The SNPs were genotyped using KASPar genotyping 
chemistry (KBioscience Ltd., Hoddesdon, UK) using the 
Fluidigm (Fluidigm Gorp., South San Francisco, GA) nano- 
fluidic 96.96 dynamic array [48], following the methods 
described by Maughan et al. [49]. End-point fluorescent 
images of the Fluidigm chips were obtained on an EP-1 
imager (Fluidigm Gorp., South San Francisco, GA). The 
data were analyzed using Fluidigm SNP genotyping ana- 
lysis software [50] . 



SNP diversity data analysis 

To visualize the ability of the SNPs to discern among cut- 
throat trout lineages, three analyses were performed. First, 
a NeighborNet phylogenetic network was created using the 
software program SplitsTree v.4.12.6 [51]. Second, principal 
coordinates analysis (PGoA) was performed using the 
Ecodist package v. 1.2.7 [52] available in the statistical soft- 
ware program R v. 2.14.1 [53]. Third, the population genetic 
software program Structure v.2.2.3 [54] and Structure 
Harvester [55] were used to infer the number of distinct 
populations within our panel. Structure analysis was evalu- 
ated 30 times for each K (with K ranging from 2 to 18), 
with 1,000,000 repetitions per run after discarding an initial 
100,000 repetitions as burn-in. In three instances, Structure 
was unable to differentiate among some subspecies (see 
Results). Each of those three groups of subspecies were 
extracted from the original input file and re-evaluated 20 
times for each K (with K ranging from 1 to 6) following the 
same procedure as outiined above. 
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MtDNA Re-sequencing 

The mitochondrial ND2 gene was re-sequenced for some 
individual samples that did not group with their pre- 
defined subspecies as expected (see Results). Amplification 
of ND2 was achieved using PCR primers Gln56F (5'-ACT 
ACA CCA CTT TCT ACT AAG GTC AGC-3') and 
AlalSR (5'-GCA TTC AGA AGA TGT GGG ATA AAG 
TC-3'). Reaction cocktails were 12.5 |il in volume, and 
contained -100 ng genomic DNA, 2.25 [il nuclease free 
water, 0.5 |il each primer, and 6.25 |^1 Promega GoTaq" 
Hot Start Green Master Mix. The thermal profOe con- 
tained an initial denature of 95°C for two minutes to acti- 
vate the enzyme, 35 cycles of 95°C for 30 seconds, 48°C 
for 30 seconds, and 72°C for 90 seconds, followed by a 
rapid cool down to 12°C. The light and heavy strands 
were each sequenced in 10.5 |il reactions using the same 
primers that amplified the ND2 gene and Big Dye chemis- 
try. Excess dye terminator was removed using Sephadex 
columns. Sequencing was performed on an ABI 3730x1 
automated sequencer located in the DNA Sequencing 
Center at Brigham Young University. 

Results 

454-Pyrosequencing, genome assembly and SNP 
discovery 

454-pyrosequencing produced 1,499,670 reads, with an 
average read length of 379 bp, for a total of 569,060,777 bp 
from a single run. The reads have been made publicly 
available via the NCBI Sequence Read Archive (Study 
#SRA062178). Reads were not equal for each subspecies. 
This inequality was likely influenced by the differences in 
sample size for each subspecies since N ranged from 2 to 4 
for each subspecies of cutthroat that was included on the 
454 plate (Table 1). Additionally, some MID-barcodes pro- 
duced more reads than others (Figure 1) even though an 
attempt was made to mix the samples in equimolar 
amounts before sequencing. The discrepancy is likely due 
to difficulties associated with fluorometric quantification of 
the PCR samples before pooling and/or inaccurate pipet- 
ting during the pooling process. Of the cutthroat trout 
lineages. Bear River cutthroat trout had the most reads 
(187,976) and coastal cutthroat trout had the fewest 
(14,607) (see Figure 2). 

We discovered 43,558 putative SNPs, of which 6383 
(15%) had BLAST hits to the GenBank refseq_protein 
database, and 37,175 (85%) did not. Our selection process 
revealed that 28,887 of those putative SNPs met the first 
five criteria oudined above for designing SNP primers 
(unique to a single subspecies, >8X coverage at the SNP, 50 
bp flanking the SNPs for primer binding sites, minor allele 
frequency of at least 3 reads, and 95% identical for a sub- 
species). Pairwise comparisons showing the number of pu- 
tative SNPs between each subspecies are shown in Table 2. 
A total of 8,627 of these putative SNPs that met the first 



five criteria were excluded from primer picking because 
they were determined to occur in repeat regions or in the 
mitochondrial genome. Of the remaining 20,260 putative 
SNPs, we developed primers for 288 SNP loci, taking care 
not to choose primers from the same contig when possible 
in an attempt to avoid picking linked loci. 

SNP genotyping and diversity analysis 

A total of 125 of the subset of 288 SNP loci for which we 
developed primers (43%) produced clean amplification 
signal, yielding genotypic clusters that were visibly sepa- 
rated from each other and could be scored following 
PCR amplification and SNP genotyping. Sequence infor- 
mation for all 125 validated SNPs has been deposited in 
the GenBank database [GenBank accession #s are listed 
in Additional file 1]. All 125 validated SNPs are listed 
with the minor allele frequencies and the subspecies that 
carried the minor allele in Additional file 2. Westslope 
cutthroat trout had the lowest number of polymorphic 
SNPs with 20, whereas greenback cutthroat trout had the 
highest number of polymorphic SNPs with 75 (Table 3). 
Westslope cutthroat trout also had the lowest number of 
highly polymorphic SNPs with 7, and Bear River cut- 
throat had the highest number of highly polymorphic 
SNPs with 22 (Table 3). 

The diversity panel was comprised of 95 individuals 
representing nine lineages of cutthroat trout (n=10 per 
lineage), rainbow trout (n=5), and a negative control. The 
total number of polymorphic SNPs ranged from 19 to 75 
for these cutthroat trout subspecies and rainbow trout 
(Table 3). Minor-allele frequencies (MAF) ranged from 
0.05 to 0.50 for each subspecies of cutthroat trout, and 
from 0.10 to 0.50 for rainbow trout (Table 3). The mean 
MAF value for all subspecies of cutthroat trout was 0.198 
per SNP locus. Because SNPs are biallelic markers, the 
maximum MAF for any given SNP locus is 0.50, which 
occurs when both SNP alleles are present at equal frequen- 
cies in the sample population. Therefore, considering SNP 
loci that exhibited a MAF > 0.3 to be highly polymorphic, 
the number of highly polymorphic SNPs in these cutthroat 
trout subspecies ranged from 7 to 22 (Table 3). Minor al- 
lele frequencies for all 125 SNPs reported here are listed in 
Additional file 2. 

The SplitsTree analysis produced a NeighborNet 
phylogenetic network that separates the cutthroat trout 
subspecies (Figure 3). Coastal cutthroat trout clustered 
together with rainbow trout. Westslope cutthroat clus- 
tered with the coastal-rainbow group, but with a large 
genetic distance between them. Lahontan Basin cut- 
throat trout exhibited great genetic distance between 
them and the other subspecies, as did the coastal- 
westslope-rainbow group. Of the interior cutthroat trout 
subspecies (i.e.. Bear River, Bonneville, Colorado River, 
Greenback, Rio Grande and Yellowstone), Bonneville 
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Figure 1 Individual read numbers. Number of reads for each cutthroat trout and rainbow trout individual, identified by BYU # and 
MID-barcode, BYU #s are prefaced with a corresponding subspecies abbreviation as follows: BR=Bear River cutthroat, BON=Bonneville cutthroat, 
COA=Coastal cutthroat, CR=Colorado River cutthroat, GB=Greenbacl< cutthroat, LAH=Lahontan Basin cutthroat, RG=Rio Grande cutthroat, 
WS=westslope cutthroat, YS=Yellow5tone cutthroat, and RBT=rainbow trout. 



and Rio Grande cutthroat separated from the others, al- 
though Rio Grande cutthroat appear in two separate 
parts of the network (Figure 3). Bear River and Yellow- 
stone cutthroat cluster together, as do Colorado River 
and greenback cutthroat trout. Moreover, 6 of the 95 
individuals did not separate as expected according to 
their a priori subspecies designation. One Lahontan 
Basin cutthroat individual (LAH95362) grouped with 



rainbow trout rather than with the other Lahontan Basin 
individuals. Three Rio Grande cutthroat individuals 
(RG90712, RG90714, and RG90732) clustered with the 
group containing Yellowstone and Bear River cutthroat ra- 
ther than with the other seven Rio Grande cutthroat indi- 
viduals. One Bear River cutthroat individual (BR239318) 
grouped with Bonneville cutthroat, and one Colorado 
River cutthroat individual (CR134518) connected to the 
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Figure 2 Read numbers by subspecies. Pooled number of reads per lineage of cutthroat trout and rainbow trout. 
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Table 2 Numbers of putative SNPs 





Bear River 


Bonneville 


Coastal 


Colorado River 


Greenback 


Lahontan Basin 


Rio Grande 


Westslope 


Yellowstone 


Bear River 


- 


2279 


482 


2170 


2430 


5459 


1807 


3600 


941 


Bonneville 


280 


- 


393 


1182 


1368 


3528 


1118 


2268 


1476 


Coastal 


50 


63 




314 


479 


407 


242 


228 


302 


Colorado River 


270 


148 


35 




1127 


2775 


969 


1679 


1225 


Greenback 


294 


131 


94 


99 




3604 


1125 


2578 


1603 


Lahontan Basin 


854 


428 


74 


358 


499 




2258 


2722 


3277 


Rio Grande 


273 


91 


24 


64 


63 


235 




1401 


1141 


Westslope 


539 


243 


50 


274 


227 


301 


152 




1930 


Yellowstone 


110 


129 


56 


110 


197 


419 


130 


200 





Pairwise comparisons of putative SNP loci among cutthroat trout subspecies discovered by 454 pyro-sequencing, Comparisons are shown for 6X coverage 
(above diagonal), and for 10X coverage (below diagonal). 



base of a branch leading to a group containing seven Rio 
Grande cutthroat individuals. Re-sequencing of the mito- 
chondrial gene ND2 showed that the "problematic" 
Lahontan Basin cutthroat individual carried rainbow trout 
mtDNA, and that the Colorado River cutthroat individual 
carried Greenback cutthroat trout mtDNA. The Rio 
Grande and Bear River cutthroat individuals carried 
mtDNA sequences of their respective subspecies. 

Principal coordinates analysis results showed that Princi- 
pal Coordinate 1 explained 35.8% of the total variance, and 
Principal Coordinate 2 explained 15.7% (Figure 4). Thus, 
the first two principal coordinates combined explained 
51.5% of the total variance observed in the distance matrix. 
Groups on the PCoA plot (Figure 4) replicate the groups 
observed on the phylogenetic network (Figure 3). 

Structure analysis and the results of Structure Harvester 
using the Evanno method revealed six distinct populations 
of cutthroat trout and rainbow trout, not the ten that we 
defined a priori (see Figure 5). Structure grouped Colorado 
River cutthroat trout and greenback cutthroat trout to- 
gether rather than as distinct subspecies, and did the same 
for Bear River and Yellowstone cutthroat. Similarly, coastal 
cutthroat, westslope cutthroat and rainbow trout were also 
considered to be a single population in the Structure 
results. Several individuals, including those that did not 
cluster as expected on the phylogenetic network, showed 



evidence of hybridization, as illustrated by mixed colors in 
the bars representing those individuals (Figure 5). Further 
Structure analyses revealed two distinct populations within 
the westslope/coastal/rainbow group, with westslope cut- 
throat cleanly separating from coastal cutthroat and rain- 
bow trout (Figure 6). Reanalysis of the Bear River/ 
Yellowstone group resulted in two distinct populations of 
cutthroat trout, although the boundary between them was 
not as clean cut as the boundary between westslope cut- 
throat and the coastal cutthroat/rainbow trout group 
(Figure 6). Reanalysis of the Colorado River/Greenback 
group resulted in three distinct populations, although many 
individuals showed signs of genetic admixture, and all indi- 
viduals that fell into the third group showed signs of 
hybridization with other subspecies of cutthroat trout 
(Figure 5 and Figure 6). 

Discussion 

We have implemented genomic and bioinformatic proto- 
cols to discover over 28,000 putative SNPs among cutthroat 
trout subspecies. We were able to scrutinize these data to 
develop a SNP assay that contains 125 nuclear SNPs and is 
capable of differentiating most subspecies of cutthroat trout 
from one another and from rainbow trout. The SNP assay 
is a fast and cost-effective way to identify individuals of un- 
known genetic background to subspecies, and can be a 



Table 3 SNP genotyping and diversity results 





Bear 
River 


Bonneville 


Coastal 


Colorado 
River 


Greenback 


Lahontan 
Basin 


Rio 

Grande 


Westslope 


Yellowstone 


Rainbow 


Sample size 


10 


10 


10 


10 


10 


10 


10 


10 


10 


5 


Polymorphic SNPs 


67 


55 


41 


65 


75 


45 


57 


20 


40 


19 


MAP range 


0.05 - 
0.50 


0.05 - 0.50 


0.05 - 
0.50 


0.05 - 0.50 


0.05 - 0.50 


0.05 - 0.50 


0.05 - 
0.50 


0.05 - 0.50 


0.05 - 0.50 


0.10- 
0.50 


Average MAP 


0.125 


0.080 


0.074 


0.130 


0.116 


0.071 


0.099 


0.045 


0.091 


0.048 


Highly polymorphic 
SNPs 


22 


11 


11 


20 


17 


11 


14 


7 


19 


10 



Summary of the results from the single nucleotide polymorphism (SNP) genotyping and diversity analysis Listed are sample size, total number of polymorphic 
SNP loci, range of minor-allele frequencies (MAP), average MAP, and total highly polymorphic SNP loci. 
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RBT94217 

LAH95362 
RBT93748 
RBT1 42374 " 
RBT1 87968 
RBT98251 - 
COA91953 COA90938 
COA92742 
COA1 30929 - 

COA91991 C°^9J504 
COA90946 
COA 130906 
■vCOAl 30905 



L AH 904 14a 
LAH90414b 



LAH 182066 
LAH91109 , 

LAH91162 
LAH90741 LAH97748 
LAH96517 
LAH96518 




BON38807 BON239318 
BON93975 BON 181 774 

BON181736 BON134577 ' 
BON 179753 



BR1 36684 BR1 79602 
YS179144 
YS1 33863 

BR1 36733 
YS141946 YS179112 
YS90212 
YS136834> 
YS1 36668 

YS136783> 

YS1 42366 



RG90714 



RG90732 
RG90712 

Figure 3 Phylogenetic network. NeighborNet phylogenetic network illustrating genetic distances among cutthroat trout species and rainbow 
trout, Individual names are represented by abbreviations for species/subspecies designations followed by BYU ID numbers. Subspecies 
abbreviations are as follows: BR=Bear River cutthroat, BON=Bonneville cutthroat, COA=Coastal cutthroat, CR=Colorado River cutthroat, 
GB=Greenback cutthroat, LAH=Lahontan Basin cutthroat, RG=Rio Grande cutthroat, WS=westslope cutthroat, YS=Yellowstone cutthroat, and 
RBT=rainbow trout. Individuals are highlighted with the same colors used to designate unique groups in the Structure results 
(see Figures 5 and 6). 



valuable tool for management agencies in their efforts to 
evaluate the genetic structure of cutthroat trout populations 
in western North America, especially prior to constructing 
and implementing conservation plans. 

These 125 putatively unlinked nuclear SNPs also allow 
for the detection of hybrid individuals, as evidenced by the 
results of our SNP diversity analyses. Indeed, the majority 
of the individuals that did not cluster as predicted in the 
phylogenetic network (see Results and Figure 3) carried a 
signature of genetic admixture between cutthroat trout 
subspecies in the results of Structure analysis (see Figure 5). 
For example, the three Rio Grande cutthroat individuals 
that clustered with Bear River and Yellowstone cutthroat in 



the phylogenetic network (RG90712, RG90714, and 
RG90732) carried Rio Grande cutthroat mtDNA haplo- 
types, but appear to carry both Rio Grande and Bear 
River- Yellowstone cutthroat alleles in the Structure bar 
graph (Figure 5). Moreover, the Lahontan Basin cutthroat 
trout individual that clustered with rainbow trout on the 
network did, in fact, carry rainbow trout mtDNA, and has 
mostly rainbow trout SNP alleles (Figure 5), thus illustrat- 
ing a case where a probable misidentified specimen was 
successfully identified using the SNP markers (the sample 
was provided as a fin clip, so we were unable to go back to 
the voucher specimen to reassess the species identification). 
The Structure analysis was also useful in detecting other 
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Figure 4 Principal coordinates analysis plot. Results of PCoA showing the first two principal coordinates that, combined, explain 51.5% of the 
observed genetic variation. Shapes on the figure legend correspond to the a priori subspecies designations. Subspecies are abbreviated the same 
as in Figures 1 and 3. Colors correspond to the population assignments made by Structure analyses (see Figures 5 and 6). 



heterozygous/introgressed individuals, even if they did 
cluster with the pre-assigned subspecies on the network, as 
evidenced by a number of mixed bars representing Bear 
River, Bonneville, coastal, Colorado River, greenback, 
Lahontan Basin and Rio Grande cutthroat trout individuals 
(Figure 5). 

The groups on the phylogenetic network (Figure 3) and 
in the PCoA results (Figure 4) accurately reflect what is 



known about phylogenetic relationships among cutthroat 
trout subspecies. It is generally accepted that coastal cut- 
throat trout was the first to branch off from the other cut- 
throat trout after the initial divergence between cutthroat 
and rainbow trout [25], so the apparent close genetic dis- 
tance between coastal cutthroat trout and rainbow trout is 
not surprising. While there are many SNPs that differenti- 
ate cutthroat trout from rainbow trout [28,40,41], the 
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Figure 5 Structure results. Results of Structure analysis showing six unique populations of cutthroat trout and rainbow trout, Groups are shown 
as follows, beginning at the upper left and moving to the right: ColoradoRiver/greenback cutthroat (red), westslope/coastal/rainbow trout (green). 
Bear River/Yellowstone cutthroat (blue), Lahontan Basin cutthroat (yellow), Bonneville cutthroat (purple), and Rio Grande cutthroat (orange). Bars 
represent unique individuals, each of which is labeled using an abbreviation for the a priori subspecies designation (given in Figure 3) followed 
by that individual's BYU identification number. 



Houston et ot. BMC Genomics 2012, 13:724 
http://www.biomedcentral.com/1471-2164/13/724 



Page 13 of 16 



— r- 



ro ^ LO (Ti 

rsl (N (N 

0 O O CTJ 

01 O CTl 

m ro fo 

rsi rsj rsi 



^iri'-i>Orviro^rNi(TiLnoo^cO'X>i^Ln:r 

(N(NOL0O'*'^Ol£>r^^O'^O 

■ — o <~i rN m 



» -J- 5 " ' 1 N ui ' : wi wi ^-f ' 7- ^Ji ut Z-i 



CT, ai in 2 
— O ^ 



^ < 
o o 



to CO g 



8 ^ 



o R 



fN (N 
00 

I- (- 

CQ CQ 



l£)rsloDoO!3:rorN^^r\lv£)fr)OOrNj^roOrnr\]iD 
(ji in oi CT) ai CTi ID — — ~ 

cn 1^ ca cn >- 



ID 



m CO 



CTi IX) wj ;^ yj ii» 

1^ 1X1 m u3 ^ cn <N 

^ \- m m m ^ 

>-'">- cQ>->-CDcQ> 



d ^ 




I I I I I I I I [ I I I I 
{Nroinoomoo*x>rvjooi^Omr\it^ 
<NO'X)mt'mrvjcji^rnr-ooooo 



in LD 



ID in (N 00 00 

00 rsi o 00 

■ ^ O O CTl ID 

^ ai o e; ^ 



O 15 ^ 



ce: q: Qc: u_ u_ ,— ^ — — 



ti_ II u_ nr Lt. 

O ^ 13 5 O U 



Figure 6 Additional Structure results. Results of secondary Structure analyses showing additional substructure in three groups that clustered as 
unique populations in the first set of analyses. Those groups are as follows: Westslope/Coastal/Rainbow, Bear River/Yellowstone, and Colorado 
River/Greenback. Secondary Structure analyses show two distinct populations in the westslope/coastal/rainbow group, two distinct populations in 
the Bear River/Yellowstone group, and three distinct populations in the Colorado River/greenback group. 



majority of the SNP primers examined herein were chosen 
specifically to detect differences among subspecies of cut- 
throat trout. Of the 288 SNP loci for which primers were 
developed in this study, only six were selected that should 
have detected differences between rainbow and cutthroat 
trout, and not aU of those amplified reliably (Additional file 
2), so the inability to clearly differentiate the two species 
likely results from an ascertainment bias that is a direct 
result of that under-sampling. 

It is somewhat surprising that the initial Structure ana- 
lysis was unable to separate westslope cutthroat from 
coastal cutthroat and rainbow trout because both the 
phylogenetic network (Figure 3) and the PCoA results 
(Figure 4) show what appears to be a large genetic distance 
between westslope cutthroat and the other lineages. The 
fact that westslope cutthroat did separate from coastal cut- 
throat and rainbow trout in the second Structure analysis 
suggests that this might have been caused by the signal 
from a small number of alleles that were unique to 



westslope cutthroat trout (Additional file 2) being overrid- 
den by the signal from the larger data-set. Lahontan Basin 
cutthroat trout also exhibit large genetic distances between 
them and other cutthroat trout subspecies on our network 
(Figure 3), and they separated from the other subspecies in 
PCoA (Figure 4) and Structure analyses (Figure 5), which is 
consistent with previously published genetic distances and 
their hypothesized positions in published phylogenies 
[20,21,25,32]. The relatively smaller genetic distances be- 
tween the other five lineages of cutthroat trout is also con- 
sistent with previously published data, and seems to 
correspond well with the evolution of these lineages in sep- 
arate drainage basins. However, results of our Structure 
analyses show that these SNPs were initially unable to dif- 
ferentiate Bear River cutthroat from Yellowstone cut- 
throat, and only did so when the Bear River/Yellowstone 
subset was reanalyzed. The inability of even the second 
Structure analyses to cleanly separate these two subspecies 
may be because the final separation between Bear River 
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and Yellowstone cutthroat likely corresponds with late 
Pleistocene events that resulted in the capture of the Bear 
River into the Bonneville Basin [26]. Gene flow was likely 
possible between Bear River and Yellowstone cutthroat 
trout populations up until the Bear River was diverted into 
the Bonneville Basin in the late Pleistocene [33-35]. It is 
possible that these two lineages have not been separated 
long enough for mutations to become fixed in each sub- 
species. It is also possible that our results are confounded 
by widespread stocking of Yellowstone cutthroat trout by 
management agencies prior to the recognition of unique 
lineages. Unfortunately, we are unable to distinguish be- 
tween these two scenarios. Similarly, we were not able to 
cleanly differentiate Colorado River and greenback cut- 
throat trout from each other using this suite of SNPs, even 
after reanalyzing the reduced datasets using Structure. It 
is unclear whether the inability of these SNPs to differenti- 
ate between Colorado River and greenback cutthroat is 
because these subspecies have diverged too recently so 
alleles have not had time to reach fixation, or if there have 
been introductions resulting in introgressive hybridization 
in what we initially treated as non-admixed populations. 
Considering the close proximity of the drainages in which 
these subspecies reside and the frequency at which cut- 
throat trout were stocked in the past, the latter scenario is 
certainly plausible. A number of bars on the Structure bar 
graph represent Colorado River and greenback cutthroat 
trout individuals that appear to be hybrids (Figure 5), which 
lends support to the latter scenario. Clearly additional re- 
search is needed to resolve this issue. 

Additional studies focused on SNP development in the 
Interior group of cutthroat trout (i.e., Bear River, Bonneville, 
Colorado River, Greenback, Rio Grande, and Yellowstone) 
are warranted, particularly when it comes to searching for 
frxed alleles between Bear River and Yellowstone cutthroat, 
and between Colorado River and greenback cutthroat trout, 
if they exist. 

Conclusions 

The SNP markers reported here have added to a rapidly 
growing body of markers that can be used in cutthroat 
trout population genetic studies, and should be a valuable 
resource in future attempts to evaluate the genetic compos- 
ition of cutthroat trout populations in western North 
America, including the detection of hybrids. These results 
reiterate that cutthroat trout subspecies are geographically 
and evolutionarUy distinct, and ought to continue to be 
managed as such by state and federal agencies. 

The method used to discover these SNP loci in cutthroat 
trout was developed for SNP discovery in the Eudicot genus 
Amaranthus [6]. Because the method was also successful 
for SNP discovery in something as evolutionary distant as 
cutthroat trout, it should be applicable to SNP discovery in 
many different kinds of non-model organisms. 



Additional files 



Additional file 1: SNP primer table, SNP marker names, GenBank 
accession numbers, the type of polymorphism for each SNP, allele 
specific primers, common reverse primers and specific target 
amplification primers are listed herein. 

Additional file 2: Characterization of each SNP locus. Major and 
minor alleles for all 125 SNP loci, along with the proportions of 
individuals within each a priori designated subspecies that carry the 
minor allele, as well as minor allele frequencies for each SNP locus 
are listed herein. 
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